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ABSTRACT 

We report the discovery of a new hysteresis effect in black hole X-ray binary state 
transitions, that of the near-infrared (NIR) flux (which most likely originates in the 
jets) versus X-ray flux. We find, looking at existing data sets, that the infrared emis- 
sion of black hole X-ray transients appears to be weaker in the low/hard state rise 
of an outburst than the low/hard state decline of an outburst at a given X-ray lumi- 
nosity. We discuss how this effect may be caused by a shift in the radiative efficiency 
of the inflowing or outflowing matter, or variations in the disc viscosity or the spec- 
trum/power of the jet. In addition we show that there is a correlation (in slope but 
not in normalisation) between infrared and X-ray luminosities on the rise and decline, 
for all three low-mass black hole X-ray binaries with well-sampled infrared and X- 
ray coverage: Lnir oc L^^^°'^. In the high/soft state this slope is much shallower; 
-^NiR (X L^^^^'^, and we find that the NIR emission in this state is most likely domi- 
nated by the viscously heated (as opposed to X-ray heated) accretion disc in all three 
sources. 

Key words: accretion, accretion discs, black hole physics. X-rays: binaries 
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1 INTRODUCTION 

X-ray binary outbursts occur when there is a substantial in- 
crease of matter accreted from the outer accretion disc or 
companion star towards the compact object (a black hole 
or a neutron star). Spectral states in X-ray binaries were 
first discovered to exist in the early 1970s (Tananbaum et 
al. 1972), when Cygnus X-1 showed a sharp transition in its 
spectral energy distribution from a state where its luminos- 
ity was dominated by soft 7— ray emission to one where it 
was dominated by soft X-ray emission. It was also found that 
this transition corresponded to a turn-off in the observed ra- 
dio emission from Cygnus X-1 (Tananbaum et al. 1972), and 
it has since been shown that this strong correlation between 
X-ray spectral shape and strength of radio emission is ubiq- 
uitous (e.g. Fender et al. 1999). Furthermore, the shape of 
the X-ray spectrum of a source is strongly correlated with 
its timing properties, with strong variability generally cor- 
related with strong radio and hard X-ray emission (see e.g. 
van der Klis 2006 for a review). 

It was observed from early on that the different spec- 
tral states occur at different fractions of the Eddington lu- 
minosity (see e.g. Nowak 1995 for a review). The low/hard 



state, which exhibits strong radio emission, occurs at sys- 
tematically lower luminosities than the high/soft states 
with low amplitude variability and no radio emission. Ob- 
servations which fit neither of these states show simi- 
lar phenomenology over a rather wide range of luminosi- 
ties; they have typically been called intermediate states 
when at lower luminosities and very high states when at 
high luminosities, but their spectral and variability prop- 
erties a re sometimes comparable and sometimes varied 
(see e.g.lHoman et al.ll200ll : iMcClintock fc Remillardll2006l : 
iHoman fc Belloni " 20051 ) . More recently, a picture has be- 
gun to develop in which these spectral state transitions 
are hysteretic, with the transition from hard state to soft 
state occurring at higher luminosities tha n the transition 
from soft state to hard sta t e (se e e.g. [M iyamot o et al.l 
19951: iNowak. Wilms fc Dovd I2OO2I: iBarret fc Olivel I2OO2I : 



Smith. Heindl fc SwankI 12003 : iMaccarone fc Coppil |2003| ). 

At sufficiently low luminosities, it appears that only 
low/hard states can exist, but the brightest low/hard states 
are just as bright as the brightest high/soft states (Homan & 
Belloni 2005). At the same time, the spread in the different 
soft-to-hard state transition levels for a variety of sources 
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appears to be considerably smaller than the spread in the 
hard-to-soft state transition levels, with a clustering in soft- 
to-hard state transition levels at about 2% of the Edding- 
ton limit (Maccarone 2003), although there do appear to be 
occasional fainter transitions (e.g. Xue, Wu & Cui 2006); 
one source, V4641 Sgr, shows behaviour that is not well 
fit by any of t he canonical states at a luminosity of about 
2 X 10"^LEdd (|Miller et al.l '2002Y For the most recent ef- 
fort at putting together a unified picture of spectral states 
in X-ray binaries, we refer the reader to Fender, Belloni & 
Gallo (2004). More recently, it has been suggested that the 
existence of a well-established relationship between X-ray 
luminosity and radio luminosity in the low/hard state with 
-/^radio oc L^^ (Corbel et al. 2000; Gallo, Fender & Poo- 
ley 2003) requires that the low/hard state be a radiatively 
inefficient fiow, with Lx oc rn? , and with the radio emis- 
sion providing a good tracer of rh , the mass accretion rate 
l|Kording. Fender fc Migliarill2006D . 

Both the rising low/hard state and the falling low/hard 
state connect to a radiatively efficient flow in the soft state, 
where Lx oc m (Shaku ra fc Sunvaev|[l973h . albeit at differ- 
ent luminosities. Changes in the bolometric luminosity dur- 
ing the transitions between the hard and soft state are gen- 
erally less than 50 percent (see discussion in iKording et al] 
I2OO6I . altho ugh a jump of more t han 50 percent may occur in 
some cases; iKalemci et"aLll2003l '). suggesting that the radia- 
tive efficiency of the accretion flow does not change sharply 
during these transitions. The radius of the inner accretion 
disc may also not change dramatically during transitions 
(e.g. lEVontera et al. I l2001al : iMiUer et al.ll2006l : iRvkoff et al.l 
I2OO7I ) although in most cases is likely to be truncated 



Esin et al.l 2001 ; Frontera et al.l 


2001bl. 2003 


; Chatv et al.l 


20031: Yuan, Cui fc Naravan 2005 


; Yuan et al. 


I2OO7). 



The X-ray luminosities in the rising low/hard state 
and the falling low/hard state have the same dependence 
on rh, but they also join to the same radiatively efficient 
soft state (in which the X-ray luminosity has a different 
rh dependence). This could indicate different radiative 
efficiency factors in the rising and falling low/hard state 
tracks. One might then also expect parallel tracks in the jet 
power versus X-ray luminosity in the two low/hard state 
tracks, related to such efficiency differences. If some or all 
of the differences between the two tracks are related to 
differences in a, the dimensionless viscosity parameter, then 
one might expect differences between the tracks due to the 
theoretically expected a dependence in the relation between 
jet power and mass accretion rate (see Meier 2001). This 
would require a to change during the soft /intermediate 
states between outburst rise and decline. Variations in a 
during transient outbursts are expected in the context of 
the disc instability model for X-ray transient outbursts, 
and are thought to be cause d by a change in the ion- 
isatio n state of hydrogen (jMever fc Mever-Hofmeisteli 
19811: ICannizzo. Ghosh fc Wheeleij 1 19821 : 



Dubus. Hameurv fc Lasotal 2001 ). or the dominant species 



in the case of accretion discs in ultracompact X-ray binaries 
which are fed by hydrogen poor white dwarfs (Menou, 
Perna fc Hernquist 2002). 

One tentative report of parallel tracks in X-rays ver- 
sus radio does exist (Nowak et al. 2005), but in general, 
the inability of radio telescopes to respond quickly to X- 



ray transients has led to most of the radio observations of 
X-ray transients being made in the decaying parts of the 
outbursts. While this would, in principle, be the case for 
optical telescopes as well, there have been dedicated moni- 
toring campaigns with small telescopes that have, in many 
cases, detected the rises of outbursts from X-ray binaries 
in the optical and infrared before they were detected in the 
X-rays, simply by monitoring known X-ray transients in qui- 
escence and waiting for them to go into outburst (see e.g. 
Jain 2001; Jain et al. 2001). 

There is mounting evidence for the spectrum of the 
jet to extend from the radio to infrared and even opti- 
cal regimes in the lo w/hard state (and hence the rise and 
decay) from spectral (|Han fc Hiellmiiig 19921: Corbel et ahl 



200ll; Ichatv et all 120031 ) and timing l|Kanbach et al.ll20o'll 



Uemura et al.ll2004l ) studies. The near-infrared (NIR) emis- 



sion has been found to be consistent wit h optic ally thin syn - 
chrotron emission (IHoman et al. I I2OO5I : iHvnes et al.l I2OO6D 
which appears to dominate the N IR regime (but proba- 
bly not the optical) in this state llCorbel fc Fender] |2003 : 
iBuxton fc Bailvnl |2004| : iRussell et al.1 120061 ) because the 
emission drops dramatically in the soft state, when the jet is 
expected to be quenched. We can therefore use the infrared 
flux to test whether there are parallel tracks in jet power 
versus X-ray power in X-ray binaries in the hard state. 



2 THE DATA 

We searched the literature for BHXBs with well sampled 
quasi-simultaneous (At < 1 day) NIR-X-ray or radio-X- 
ray data on the rise and/or decline of a transient outburst 
(where NIR refers to J, H or K-hand). Data from three 
sources were found: 4U 1543-47, XTE J1550-564 and GX 
339-4. XTE J1550-564 is the only source for which we ob- 
tain both rise and decay data of the same outburst. Appar- 
ent fluxes were converted to intrinsic luminosities adopting 
the best known estimates of the distance to the sources, 
the extinction Ay and neutral hydrogen column density 
A^'h. In Table 1 we list these values and the outbursts, tele- 
scopes used and references of the data. We adopt Ay = 5 
for XTE J 1550-564 as oppose d to Ay = 2.5, which is 
adopted in iRussell et al. I (|2006h . A value of Ay ~ 5 is in- 
ferred from the hydrogen column density measu red from X- 
ray observations of this source tTomsick. Corbel fc Kaaretl 
I2OOII : iTomsick et all l2003l : iKaaret et al]|2003l ). A value of 
Ay ~ 2.5 pr oduces a red optical spectrum (7 < where 
oc i^^; Rus sell et al.|[2006 ). which would indicate a much 
cooler disc than is expected, whereas with Ay ~ 5 the spec- 
trum is blue (7 > 0), as is expected for the Rayleigh- Jeans 
side of the disc blackbody spectrum. While uncertainties in 
the extinction will affect the relative normalisations of the 
Lnir-Lx relations for different sources, they will not affect 
the slopes of the relations for the individual sources. 

To calculate the X-ray unabsorbed 2-10 keV luminos- 
ity we assume a power-law with a spectral index = —0.6 
(photon index F = 1.6) when the source is in the hard state 
and a blackbody at a temperature of 1 keV whe n the source 
is in the soft state (the se are the same as in ICallo et al] 
l2003l : IRusscII et al.ll2006l ). The resulting luminosities do not 
change significantly if similar reasonable approximations to 
these values are adopted. RXTE ASM X-ray counts were 
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Table 1. The sources used and data collected to calculate their luminosities. 



Source Outburst Telescopes Distance Nn / Flux 

/ kpc lO'^^ cm~^ refs 



4U 1543-47 


2002; decline 


YALO 1.0 m, RXTE PCA/HEXTE/ASM 


7.5±0.5 


1.55±0.15 


4.3±0.2 


1, 2 


XTE J1550-564 


2000; rise + decline 


YALO 1.0 m, RXTE PCA/HEXTE/ASM 


5.3±2.3 


5.0 


8.7dz2.1 


3, 4 


GX 339-4 


2002; rise 


YALO 1.0 m, RXTE PCA/ASM 


O+7.0 
*-1.0 


3.9±0.5 


g+0.9 
"-1.7 


5 


GX 339-4 


1999; decline 


ATCAi, RXTE PGA 
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Australia Telescope Gompact Array. T he distances an d extinction estimat es are from the references in Table 1 of lRussell et al ] l|200d) . 
except Ay = 5 for XTE J155 0-564 llTomsick et al.ll2001. 2003; Kaa ret e t al.ll2003l. see tex t ^. Refe ren ces for the fluxes use d: (1) 
iBuxton fc Bailvnl l|2004h : f2') lKalemci et al.l l|2005E : ( S*) fjain et"aLl 1 I2OOII) ; ('4') lMiller et"all j200lh : f5') lHoman et al.l l|2005h : (6) 

iNowak et al. 
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Figure 1. The NIR and X-ray light curves of the 2000 outburst of 
XTE J1550-564 (upper panel). The lines are fit to the initial rise 
data before the break in the rise behaviour of both light curves. 
All X-ray data are from the RXTE ASM except in the hard state 
decline, which are from the RXTE PGA/HEXTE. The ratio of 
X-ray to NIR fluxes are shown in the lower panel. 



converted to flux units using the NASA tool Web-PIMMS. 
Since we use data close to state transitions, we adopt the 
dates of these transitions as defined by the analysis of the 
authors in the literature. Inaccurate background subtraction 
could severely affect the X-r ay luminosity e s timat es at low 
luminosities. For 4U 1543-47, iKalemci et al.1 (|2005l 'l have ac- 
counted for the Galactic X-ray background. The X-ray data 
we use for XTE J1550-564 is more than two orders of magni- 



tude above the quiescent level so no background subtraction 
is necessary. The X- ray background co ntamination for GX 
339-4 is uncertain (|Nowak et al.ll2002l ), but for a reason- 
able choice of background fluxes the slope of the resulting 
NIR-X-ray correlation (see Section 3) does not change. 

For NIR data near quiescence (or close to the lowest 
NIR fluxes observed from the source) we only include data 
which are > 0.5 mag above the measured quiescent flux level. 
We also subtract this quiescent level from all NIR data so 
that only the outburst flux remains (the NIR in quiescence is 
likely dominated by the companion star) . The well-sampled 
NIR quiescent level of XTE J1550-564 is brighter by ~ 0.25 
magnitudes afte r the 2000 outburst than prior to it (Fig la of 
iJain et al ] |200ll ). This enhancement could be due to thermal 
emission from the heated accretion disc, or jet emission. 

It is possible that the star or outer regions of the ac- 
cretion disc are re-emitting X-ray radiation which was ab- 
sorbed during the outburst. The donor s tar in XTE J1550- 
564 (which has spectral type G8 to K4: lOrosz et al.ll2002l ) 
occupies around 10 percent of the solid angle seen by the 
X-ray source and shoul d absorb ~ 70 percent of the radia- 
tion incident on it (Basko &: Sunvaevl ll973). Since the out- 
burst gave off ~ 3 X 10*^ ergs and the source brightened by 
about 10'''^ ergs s~^ in the infrared, it could maintain this 
brightening for 70 years if all the energy was lost in this 
infrared waveband. The brightening was present in at least 
the optical IZ-band too, but the brightening of the broad- 
band optical-infrared spectrum could still be sustained for 
more than one year, longer than the duration of the quies- 
cent li ght curve used in this paper. In addition, lOrosz et al.l 
(|2002l ) found that the quiescent optical flux level of XTE 
J1550-564 was lower in May 2001 than before and after 
the 2000 outburst. We therefore have subtracted the (not- 
dereddened) base levels oi H = 16.25 from the data on the 
hard state rise and H = 16.0 from the hard state decline 
data. We also note that long-term variations in infrared qui- 
escent levels have been observed in GRO J0422-I-32, A0620- 
00 and Aql X-1 (Reynolds, Callanan & Fihppenko 2007; C. 
Bailyn, private communication). 

NIR intrinsic (de-reddene d using the extinction Ay 
listed in Table 1, according to ICardelli. Clayton fc MathisI 
Il989h luminosities were calculated adopting the approxima- 
tion Lnir ~ vFy^NiR (we are approximating the spectral 
range of each filter to the central wavelength of its wave- 
band; the same method as lRussell et al.ll2006l ). 
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3 RESULTS 

The X-ray and NIR light curves of the 2000 outburst of 
XTE J 1550-564 are plotted in Fig. 1 (upper panel; see also 
I Jain et aLll200l[ ). In both X-ray and NIR light curves there 
is a break in the relation between luminosity and time on 
the outburst rise. One explanation for this could be that 
the point at which the luminosity versus time curve flattens 
corresponds to the start of the hard/intermediate state, and 
not to the bona fide low/hard state. However, this tradition- 
ally occurs when the X-ray spectrum begins to soften, which 
happened about six days after the X-r ay break (the firs t 
transition dashed line in Fig. 1: see also iMiller et all l200ll ) . 
We note also that the flattening in the infrared appears to 
happen a few days before the flattening in the X-rays, and 
the optic al (V-hanA) pre cedes the NIR (i/-band) also by 
~ 3 days l|jain et al.ll200ll ). One may expect these delays to 
reveal a relation between wavelength and delay time, but the 
order between the wavebands is incorrect for this to be the 
case; the optical change is followed by the infrared and then 
X-ray. T he reasons for these d elays between wavebands are 
not clear. lUemura et al.l l|2000l ) found that the X-rays lagged 
the optical rise of an outburst of XTE J1118-I-480 by about 
10 days, and interpreted it as in 'outside-in' outburst; the 
optical emission originating from the viscously-heated (as 
opposed to X-ray heated) accretion disc. In the case here of 
XTE J1550-564, we would not expect the NIR light to lag 
the optical, as is observed, if the viscous disc is the origin of 
the emission. 

We have fit exponentials in luminosity versus time for 
the initial rise data before the break to a shallower slope in 
both light curves (Fig. 1 upper panel). These are required to 
obtain the Lnir-I/x relation before the break, as there is lit- 
tle quasi-simultaneous data on the rise but clear luminosity- 
time relations exist. For the X-ray light curve we have made 
a fit to the first four RXTE ASM data points (all have > 3o- 
signiflcance) . The break appears just after these data in the 
X-rays and all four points lie very close to the fit. All H- 
band data before the gap at MJD 51634-40 were used for 
the LNiR-time fit. 

In Fig. 2 we plot the quasi-simultaneous NIR (or radio) 
and X-ray luminosities for each source. The luminosity-time 
fits to the rise of the XTE J1550-564 outburst in NIR and 
X-ray (Fig. 1) are used to infer the Lnir-I/x relation on the 
rise for this source before the data are sampled well. This fit 
has a slope Lnir cx: where j3 = 0.7. Interestingly, fitting 
the data after the aforementioned NIR and X-ray breaks 
(and before the state transition) in the XTE J1550-564 light 
curve yields the same slope: 13 = 0.7. 

For GX 339-4 we have NIR rise data and radio de- 
cline data from different outbursts (|Homan et al.l l2005l : 
iNowak et al]|2005^ . We therefore cannot compare the nor- 
malisations of the rise and decline in this source even if 
we assume the two outbursts were similar hysteretically; 
assumptions are required of the radio-to-NIR spectral in- 
dex. For 4U 1543-47 we have outburst decline data only, 
but there is an increase in the NIR fiux when the source 
enters the hard state, as is expected from jet emission 
iHoman et al.ll2005l : iRusseU et aLll2006l ). 

We see from Fig. 2 that the NIR is enhanced at a given 
X-ray luminosity on the decline compared to the rise for the 
only source with both rise and decline data from the same 
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Figure 2. X-ray luminosity versus quasi-simultaneous monochro- 
matic NIR/radio luminosity for the three sources. Crosses are 
hard state data and filled circles represent data in the soft state 
when the jet is quenched. The NIR data are J-band for 4U 1543— 
47 and iT-band for XTE J1550-564 and GX 339-4. The data on 
the decline for GX 339—4 are radio 8.64 GHz from a different 
outburst, not NIR (and is the only radio data we use). 
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Figure 3. The slope fS of the Lnir-Lx relation for the decline 
data of 4U 1543-47 and XTE J1550-564, as a function of the 
maximum Lx of the data used to infer the fit. 
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Table 2. Measured slopes of the Lnir^^'X (or I/radio^^'x) rela- 
tions. Lnir <^ L^. 



Source 




Data 


Slope P 


4U 1543-47 




soft state 


0.13±0.01 


4U 1543-47 




hard state decline 


0.5-0.8 


GX 339-4 




hard state rise 


0.55±0.01 


GX 339-4 




hard state decline^ 


0.70±0.06 


XTE J1550- 


-564 


hard state rise 


0.69±0.05 


XTE J1550- 


-564 


soft state-*- 


0.23±0.02 


XTE J1550- 


-564 


hard state decline 


0.63±0.02 



^AU soft state data are used for this fit except those above the 
apparent relation at the beginning and end of the soft state phase 
(likely to be during state transition). ^The slope here is measured 
from the five radio data points. 



outburst, XTE J1550-564 (based on the extrapolation of 
the X-ray data on the rise, as explained above). The level of 
enhancement is ~ 0.7 dex in NIR luminosity, or a factor of ~ 
5. The NIR luminosity at a given Lx (i.e. the normalisation) 
also differs between sources; this may reflect distance and/or 
interstellar absorption uncertainties. 

A striking observation from Fig. 2 is the similarity of the 
slope of the Lnir-Lx (or I/radio-ix) relation in the low/hard 
state rise and decline between sources. In Table 2 we tabu- 
late the slopes measured for each source. The slope of the re- 
lation for 4U 1543-47 and XTE J1550-564 on their declines 
depends on which hard state data are used, since there is an 
initial NIR rise after the transition to the hard state. In Fig. 
3 we show the inferred slope /3 as a function of the range 
of data used to measure the fit. We use all the data at low 
X-ray luminosity and impose a cut at a higher X-ray lumi- 
nosity: L^''". We plot 13 versus log(Lx'"') in Fig. 3. There 
appears to be a smooth relation between /3 and ix'"'i the 
slope becomes shallower when we include the data at high 
Lx because of the initial NIR rise as the source declines in 
Lx- The cut required to measure the Lnir-Lx relation once 
the two wavebands are positively correlated should be at a 
value of Lx after the NIR begins to drop. At this cut we 
have f3 - 0.5-0.8 and /? = 0.63 ± 0.02 for 4U 1543-47 and 
XTE J1550-564, respectively. From Table 2 we can see that 
generally, the slopes in the hard state (/3 = 0.5-0.7) are in 
agreement with the radio-X-ray and optical/NI R-X-ray cor- 
relations found from these and other sources ifCorbel et al.l 
l2003l : iGallo et al.|[2003l . l2006l : iRusseU et af]|20(m 

In addition, Fig. 2 shows that all the NIR and radio 
data we use are suppressed in the soft state, indicating its 
most likely origin in the low/hard state is the jet, which is 
expected to be quenched (see Section 1). In addition, the 
optical-NIR colours are red (7 < 0) in all three sources in 



the hard state, indicating optically thin emission ( Jain et al] 
I2OOII : iBuxton fc Bailvnl 120041 : iHoman et all 120051 ). The op- 



tical and NIR luminosity from the viscously-heated accre- 
tion disc should vary relatively smoothly across the state 
transition, as should X-ray reprocessing in the disc unless it 
largely depends on the hard spectral component (see Russell 
et al. 2006 for a more detailed discussion of these and other 
sources). The emission may also not be smooth if it is corre- 
lated with a and if a changes across the transition (see Sec- 
tion 4.2). Since we see a stronger soft state quenching in the 




Mass accretion rate (arbitrary units) 



Figure 4. A schematic diagram of the X-ray luminosity versus 
accretion rate. The units of the accretion rate are such that, in 
the 'normal' low/hard state, at the state transition back to the 
low/hard state, the X-ray luminosity is 2% of the Eddington lu- 
minosity and the mass accretion rate is unity. The solid line corre- 
sponds to the decaying low/hard state, and to sources like Cygnus 
X-1 which show n o state transition/luminosity hysteresis effects 
jSmith et al.ll2002l) . while the dashed line corresponds to the ris- 
ing portion of hysteretic light curves. The arrows trace a typical 
outburst's track in this diagram. The dot-dashed line shows the 
evolution of the jet kinetic power, as inferred from the observed 
hysteresis effect. It drops below the y-axis in the high soft state 
to indicate that the jet is 'quenched' in this state. 



NIR than in the optical wavebands in our sources l|jain et al.l 
l200ll : lBuxton fc Bailvnll2004lHoman et alllioosl ). the longer 
wavelengths have much larger ratios of jet to disc contribu- 
tions (this is expected since the spectral index of the jet 
spectrum is less than that of the disc). 

In Table 2 we also state the slopes of the Lnir-Lx re- 
lations in the soft state for 4U 1543-47 and XTE J1550- 
564. The theoretically expected slope if the origin of the 
NIR emi ssion is X-ray repro cessing in the disc is f3 — 0.5 
(ivan Paradiis fc McChntocklir994) . The shallower slopes ob- 
served are close to those expected if the origi n of the NIR 
emission is the viscously heated disc: /3 = 0.15 l| Russell et al.l 
2006). Since we are using a finite X-ray energy range (2-10 
keV) as opposed to the bolometric X-ray luminosity, (3 may 
be inaccurate because the fraction of the bolometric X-ray 
luminosity in this range will change with Lx in the soft state. 
However, this fraction would have to change by a factor of 
several during the ~ one dex drop in Lx for (3 to be inaccu- 
rate by a factor of ~ 3-4 (as is required if /3 ~ 0.5 for X-ray 
reprocessing). We can therefore conclude from correlation 
analysis that the viscously heated disc is likel y to dominate 
the NIR in the soft state in these two sources. iHoman et al-l 
(2005) also found this to be the dominating NIR emission 
process for GX 339-4 in the soft state (from the same data 
that we use here) since changes in the NIR preceded those 
in the X-ray by several weeks. 

It is unlikely that disc emission contaminating the in- 
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frared luminosity plays a substantial role in the observed 
hysteresis effect in XTE J1550-564 in the hard state, be- 
cause the separation between the rise and decline parallel 
tracks is a factor of about five in Lnir; the disc component 
would have to be dominating the NIR luminosity in the ris- 
ing hard state for this to be the explanation of the effect. 
It would then be quite surprising for the colours to change 
so dramatically at the state transition and would have to 
be for some other reason than the quenching of the jet. In 
the decaying hard state of XTE J1550-564, a rise is seen 
in I and V as well as in H at the time of the state transi- 
tion, indicating again that the jet power-to-disc power ratio 
is larger for the decaying hard state than for the rising hard 
state (see the discussion of the optical 'reflare' in Jain et al. 
2001). 



4 DISCUSSION 

In this Section we concentrate mainly on interpreting the 
hysteresis effect. Nearly all models for jet production sug- 
gest that the kinetic power supplied to the relativistic jet will 
be linearly proportional to the mass accretion rate (e.g. Fal- 
cke & Biermann 1995; Meier 2001). The observed monochro- 
matic luminosity of the synchrotron self-absorbed part of the 
jet spectrum will be proportional to the kinetic luminosity of 
the jet to the 1.4 power (e.g. Blandford & Konigl 1979; Fal- 
cke & Biermann 1995; Heinz & Sunyaev 2003). As a result, 
infrared and longer wavelength fluxes may be expected to be 
excellent tracers of the mass accretion rate, and would be 
much less sensitive to c hanges in the accretio n efficiency than 
at X-ray wavelengths (|Kording et al.|[2006l ). However, in at 
least one neutron star X-ray binary 4U 0614-1-09, the NIR jet 
contribution is optically thin with the break between opti- 
cally thick and o ptically thin emission lying further into the 
mid-infrared tMigliari et al.1,200 6). For B HXBs, it has been 
shown that this break lies close to the NIR l|Corbel fc Fendeil 



l2002l: iRusseh et aJl2006l) but may change with luminosity 



( Nowak et al.ll2005 

For a jet being powered from an advection dominated 
accretion flow (ADAF; Narayan & Yi 1995), it has been 
suggested that the relation between m and the jet's kinetic 
luminosity can also depend on the dimensionless viscosity 
parameter a, and on some fudge factors related to the fact 
that the standard ADAF treatment of rotational velocities 
and azimuthal magnetic field strengths are only approxima- 
tions (Meier 2001). This assumes that an ADAF produces 
the X-ray emission. Narayan & Yi (1995) did not account 
for the fact that shear in the Kerr metric will enhance az- 
imuthal magnetic fields, even for black holes rotating at a 
fraction of their maximal rate (e.g. Meier 1999). We then 
find that: 



i-^j oc ma 



(1) 



where Pj is the total power injected into the jet, m is the 
mass accretion rate through the disc, and f{g,j) is a correc- 
tion term which is a function of the black hole spin j and 
the fudge factors which are rolled into g. 

In the following subsections we use these arguments and 
discuss possible explanations for the observed hysteresis ef- 
fect. 



4.1 Hysteresis traces the radiative efficiency of 
the accretion flow? 

Let us first consider the simplest possible interpretation of 
the parallel tracks in XTE J1550-564 - that they represent 
variations in the radiative efficiency in the rising and de- 
caying low/hard states, with the mass accretion rate being 
traced perfectly by the infrared fiux. Here, we effectively 
assume that the variations in a and the fudge factors are 
unimportant, as are variations in the radiative efficiency 
of the jet itself for producing infrared emission. Since both 
the rising low/hard state and the decaying low/hard state 
are observed to follow the I/radio oc ix"'^ relation (Gallo et 
al. 2003; 2006), they are both expected to have Lx oc 
(a radiatively inefficient fiow). At the same time, they are 
both expected to join to the same high/soft state solution, 
which is believed to be a radiatively efficient (Lx oc m) ge- 
ometrically thin accretion disc (Shakura & Sunyaev 1973). 
XTE J1550-564 was said to be in an intermediate state be- 
tween the rise and decline of the 2000 outburst based on 
timing properties and a power law component in the X-ray 
spectrum. However, the black body flux was alw ays domi- 
nant over the power law fiux (|Miller et al.ll2003l ). and the 
radio/infrared jet was quenched so for purposes of this work 
we can treat it as a high/s oft state from its jet properties 
(see e.g. iFender et al.ll20M ) in the context of our following 
analysis. A schematic version of the expected track followed 
by a hysteretic source is plotted in Fig. 4. We then see that 
the scenario where rh is well traced in the hard state by the 
infrared luminosity and in the soft state by the X-ray lumi- 
nosity is immediately invalidated as this would require the 
infrared luminosity to be lower at a given Lx on the decline 
compared to the rise - the opposite to what is observed. To 
achieve what is observed, we instead require the jet (NIR) 
luminosity to evolve as the dot-dashed lines in Fig. 4. 



4.2 Hysteresis traces q? 

Next, let us consider the possibility that rh is instead sensi- 
tive to variations in a (Pj and a are linked through equation 
1). This is a natural explanation for the different state tran- 
sition luminosities in the rising and falling low/hard states; 
it is found by Esin, McClintock & Narayan (1997), for ex- 
ample, that the state transition luminosity should be 1.3a^, 
while it is found by Zdziarski (1998) that the state tran- 
sition luminosity should scale approximately as y^^^a^^^, 
where y is the Kompaneets parameter, and is defined as 
y — 4fcT/(moC^)max(r, T^), where rric is the electron rest 
mass, k is the Boltzmann constant, T is the temperature 
of the Comptonising region, and r is the optical depth of 
the Comptonising region. In either of these cases, then, one 
would qualitatively expect that increasing a would increase 
the state transition luminosity. A similar result, albeit with- 
out an analytic approximation for the size of the state transi- 
tion luminosity, is found by Meyer-Hofmeister, Liu & Meyer 
(2005), who showed that hysteresis may result from the dif- 
ferent efficiencies of hard photons and soft photons in evap- 
orating mass from the thin disc into the corona. 

Some of the other hysteresis properties of outbursts 
seem to be qualitatively consistent with the idea that a vari- 
ations play at least some role. A larger mass density at the 
onset of an outburst may lead to a higher a when the insta- 
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bility is triggered. If this is the case, then it will naturally 
follow that both the peak luminosity of the outburst, and 
the luminosity at which the hard-to-soft transition occurs, 
will be correlated with o ne another through the connectio n 
to a, as has been shown ||Yu. van der Klis fc Fendeij|2004l ). 

The radiative efficiency of the hot disc will depend on 
which LmR-Lx track one chooses, and hence on a. In gen- 
eral, the X-ray luminosity is given by Lx = ern, where e is 
the radiative efficiency. For a solution where Lx oc m^, then 
e oc m. More specifically, when the two hard state tracks 
must join onto the same radiatively efficient soft state track, 
then e oc where rricr is the mass accretion rate at the 

TTlcr 

critical state transition luminosity (which has one value for 
the rise track and a different one for the decline). Let us 
now take a general form for the a dependence of this value, 
rrici — a'^. The radiative efficiency of the hot disc will then 
scale as Q~'^. As a result, rh at a given X-ray luminosity will 
scale as a~'^^^ , since Lx oc rr? in the hard state. 

The a variations will thus not lead to parallel tracks in 
^NiR versus Lx if the state transition luminosity scales as 
. The a dependence of the jet power extraction efficiency 
will cancel out exactly with the a dependence of the X-ray 
radiative luminosity. The estimate of the state transition 
luminosity given by Esin et al. (1997) can therefore be valid 
only if the fudge factors vary hysteretically or if the radiative 
efficiency of the jet varies hysteretically. The dependence of 
the transition luminosity on a must be weaker than o? in 
order for a variations to be solely responsible for the parallel 
tracks in Lx versus Lnir- 

The estimate of the state transition luminosity put forth 
by Zdziarski (1998) does have a weaker than o? dependence 
for the state transition luminosity: ci^'\ This leads to a re- 
lation where Lx oc L-^ctO?^^ ■, a closer match to the observa- 
tions, but still problematic quantitatively. For a given Lnir, 
the X-ray luminosity is ~ 10 times higher in the rising hard 
state than the decaying one. The difference between the X- 
ray luminosity at the hard-to-soft transition compared to 
the soft-to-hard transition is a factor of about five. If the 
observations are to be interpreted in terms of variations in 
a only, then either the state transition luminosity must scale 
slightly weaker than linearly with a, or the jet power must 
scale more strongly than linearly with a (e.g. if the fudge 
factors in the Meier 2001 relation scale with a in some way). 

4.3 Hysteresis traces jet behaviour? 

The other alternative is that the spectral shape or the radia- 
tive efficiency of the jet may change hysteretically. A test of 
whether the jet power, or the jet spectrum is the most impor- 
tant factor may come from correlated high resolution timing 
in the X-rays and optical/ infrared bands. It has previously 
been shown that a negatively lagged anti-correlation of opti- 
cal emission relative to X-ray emission (Kanbach et al. 2001; 
Hynes et al. 2003) can be explained as a result of the jet tak- 
ing away the bulk of the kinetic luminosity from the accre- 
tion flow (Malzac, Merloni & Fabian 2004). If the hysteresis 
effect we see is the result of a different amount of kinetic 
power being taken away by the jet at the same accretion disc 
rh, then there should be a signature in the magnitude of the 
negatively lagged anti-correlation. We note that observing 
the spectral changes directly would be preferable, but this 
is quite difficult due to the convergence around the J-band 



of the optical components from the accr etion disc, the donor 
star, and the jet i n most systems (e.g. iHoman et al.ll2005l : 
iRussell et aLll2006l ). 

A completely independent emission mechanism does ex- 
ist for explaining the X-ray emission in X-ray binaries - 
jet emission (e.g. Markoff, Falcke & Fender 2001; Markoff, 
Nowak & Wilms 2005). In this scenario, the observed par- 
allel tracks would form a jet-jet hysteresis as opposed to 
a jet-inflow (e.g. jet-corona) hysteresis. In the pure jet- 
synchrotron model (Markoff et al. 2001), the NIR-to-X-ray 
ratio can change at a given Lnir only if either the spec- 
tral slope changes, or the break frequency where the jet be- 
comes optically thin changes (see e.g. Nowak et al. 2005 
for a discussion of attempts to test the model based on 
observed spectral shapes). In the more recent, more com- 
plicated scenarios, where the corona forms the 'base of the 
jet' (Markoff et al. 2005) there are additional parameters 
which can change this flux ratio, by changing the fraction 
of the power that comes from Compton scattering at the 
base of the jet. Given that the physics of the underlying 
accretion flow feeding the jet is considerably less developed 
in these models than for ADAF and related models, there 
are no quantitative determinations of what should be the 
state transition luminosities or how the resulting hysteresis 
in the state transitions should behave. It thus seems prema- 
ture to comment on whether the observations presented in 
this paper have any implications for or against the viability 
of jet X-ray models, but we note that we here provide an 
additional constraint for such models. 

4.4 Consequences of the interpretation 

These findings indicate, on both observational and theoreti- 
cal grounds, that there is reason to expect substantial scatter 
in relations between radio/NIR and X-ray luminosity espe- 
cially if some of the data points are taken during the rising 
parts of outbursts. At the present time, it seems that most 
of t he scatter in t he re lations of Gallo et al. (2003, 2006) 
and I Russell et al"] (|2006l ) are due to the uncertainties in the 
measurements of the masses and distances of the different 
systems used to make the correlations, ho wever two recently 
discovere d sources could be exceptions fCadoUe Bel et al.l 
2007; Ro driguez et al.ll2007l : see also [Callo 200"^). As more 
data is collected for some of the sources (especially GX 339- 
4), it may become clear that the relation is multiple- valued, 
and the hysteresis effects will need to be considered. 



5 CONCLUSIONS 

We have shown that there is an intrinsic correlation between 
NIR (which we show comes from the jets) and X-ray lumi- 
nosity; Lnir oc L^^""'^ for all three low-mass black hole 
X-ray binaries in the low/hard state which have well sam- 
pl ed data. This is not just a global correlation as shown 
in iRussell et al.l l|2006l ). but a correlation that holds in all 
hard state rise and declines. In the high/soft state, there is 
evidence for the NIR emission to be dominated by the vis- 
cously heated disc in all three sources, with Lnir oc L^^""'^ 
observed in two sources. 

The normalisation of the hard state correlation can dif- 
fer with time for a single source. The infrared flux in the 
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2000 outburst of XTE J1550-564 shows a hysteresis eflfect, 
in that the rising low/hard state has a lower infrared flux at 
a given X-ray flux than the decaying low/hard state. Several 
theoretical explanations for this effect seem plausible, and 
have direct implications for the accretion process, jet power 
and the dominating X-ray emission processes. We have ruled 
out the possibility that the hysteresis is directly caused by 
different radiative efficiencies in the inner accretion flow on 
the rise and decline. The effect may be due to changes in 
the viscosity parameter a and/or the spectrum/radiative 
efBciency of the jet. Hopefully, these results will help mo- 
tivate future observational searches for more examples of 
jet-disc hysteresis effects and more theoretical work which 
might help to understand whether the disc viscosity param- 
eter ct should vary as a function of time over X-ray binary 
outbursts. 
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